Experimental study of decays to K^K 7r° and K'^K rj 
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Abstract 

Using (106 lb 4) X 10^ V' events accumulated with the BESIII detector at the BEPCII e+e~ 
collider, we present measurements of the branching fractions for ij/ decays to K^K^ir^ and 
K+K-f]. In these final states, the decay if^' K^{U30)+K- + c.c. is observed for the first 
time, and its branching fraction is measured to be (7.12 it 0.62 (stat.)lQ'gJ (syst.)) x 10~^, which 
indicates a violation of the helicity selection rule in tp' decays. The branching fractions of ip' 
K*{892)'^K~ + CO., (jn], 4>Tr^ are also measured. The measurements are used to test the QCD 
predictions on charmonium decays. 

PACS numbers: 13.25.Gv, 13.20.Gd, 14.40.Pq 
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I. INTRODUCTION 



In the framework of perturbative QCD (pQCD), J /ip and ip' decays to light hadrons are 
expected to be dominated by the annihilation of cc quarks into three gluons or one virtual 
photon, with hadron decay partial widths that are proportional to the square of the cc wave 
function overlaps at the origin, which can be related to their leptonic decay widths [l| . This 
suggests that the ratio Qh of branching fractions for ip' and J /ip decays to the same final 
state should follow the rule: 

^'^ Br{J/^ ^h)~ Br{J/^ ^ e+e") ~ ' 

where Br denotes a branching fraction and /i is a particular hadronic final state. This 
relation is referred to as the "12% rule" . 

Although the 12% rule works well for some specific decay modes of the the decay 
ip' to pvr exhibits a factor of 70 times stronger suppression than expectations based on 
this rule. This suppression in vector-pseudoscalar (VP) meson modes was first observed by 
MARKII [2,], which is referred to as the "pvr puzzle". Further tests of this rule in the VP 
modes have been performed by CLEO [3] and BESII [4], and have been extended to the 
pseudoscalar-pseudoscalar meson (PP), vector-tensor meson (VT) and multibody decays. 
Although Qh values have been measured for a wide variety of final states, most of them 
have large uncertainties due to low statistics jsl. Reviews of the rho-pi puzzle conclude 
that current theoretical explanations are unsatisfactory 6|. More experimental results are 
desirable. 

For charmonium ip{X) decays to light hadrons hi{Xi) and h2{X2), the asymptotic behavior 
of the branching fraction from a pQCD calculation to leading twist accuracy gives [3]: 

5r[^(A) ^ /.i(A0/^2(A2)] ~ ' , (2) 

where A, Ai and A2 denote the helicities of the corresponding hadrons. Here rric is the 
charm quark mass and Aqcd is the QCD energy scale factor. If the light quark masses are 
neglected, the vector-gluon coupling conserves quark helicity and this leads to the helicity 
selection rule (HSR) [8|: Ai + A2 = 0. If the hehcity configurations do not satisfy this 
relation, the branching fraction should be suppressed. 
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For the V' decays to VP [K*{892)^K^] or TP [K^{U30)^K^], the amphtudes are anti- 
symmetric in terms of the final state hehcities, since strong or electromagnetic interactions 
conserve parity. Hence the amplitudes vanish when Ai = A2 = 0. Nonvanishing amplitudes 
require the helicity configuration to satisfy the relation |Ai + A2I = 1, which violates the 
HSR and the branching fractions are expected to be suppressed. 

Strikingly, HSR-violating decays were recently observed in Xcj decays into vector-vector 
meson pairs by BESIII jsl, which strongly indicates the failure of the HSR jioj]. In an analysis 
of ifj' KlK^-K^ by BESII fl, evidence for i)' K*jK^ {K} refers to either K*j{U30) or 
i^*(1410)) was seen, but low statistics prevented a further study. 

With the large ip' data sample accumulated by the BESIII experiment, new opportunities 
to precisely test the 12% rule in the decays of ip' — ^ K* (892)^ K~ + c.c. and rjcj), and to search 
for ip' — )■ K2{14:30)^K^ are available. Such measurements can shed light on charmonium 
decay mechanisms and, therefore, be helpful for understanding the pn puzzle. In particular, 
the decay ip' — > K^K^rj provides opportunities to study not only 0?7, but also the excited 
states, such as 03(1850) and 0(2170). The decay tp' — )■ K^K'i:^ also allows us to study the 
isospin violation decay ip' 07r°, which is expected to proceed via electromagnetic (EM) 
processes jll |. 



II. THE BESIII EXPERIMENT AND DATA SET 



We use a data sample containing (106 ± 4) x 10^ ip' decays recorded with the BESIII 



detector [l2| at the energy-symmetric double ring e^e~ collider BEPCII. The primary data 
sample corresponds to an integrated luminosity of 156.4 pb~^ collected at the peak of the %p' 
resonance. In addition, a 2.9 fb~^(43 pb^^) data sample collected at a center-of-mass energy 
of 3.773 GeV (3.65 GeV) is used for continuum background studies. 

BEPCII is designed to provide a peak luminosity of 10^'^ cm~^s~^ at a beam current of 
0.93 A for studies of hadron spectroscopy and r— charm physics [13] . The BESIII detector 



is described in detail elsewhere [12j. Charged particle momenta are measured with a small- 
celled, helium-gas-based main drift chamber (MDC) with 43 layers operating within the 
IT magnetic field of a solenoidal superconducting magnet. Charged particle identification is 
provided by measurements of the specific ionization energy loss dE / dx in the tracking device 
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and by means of a plastic scintillator time of flight (TOF) system comprised of a barrel part 
and two endcaps. Photons are detected and their energies and positions measured with 
an electromagnetic calorimeter (EMC) consisting of 6240 CsI(Tl) crystals arranged in a 
barrel and two endcaps. The return yoke of the magnet is instrumented with resistive plate 
chambers arranged in 9 (barrel) and 8 layers (endcaps) for the discrimination of muons and 
charged hadrons. 

The optimization of the event selection criteria and the estimation of background sources 
are performed with Monte Carlo (MC) simulated data samples. The GEANT4-based simula- 
tion software [l^ includes the geometric and material description of the BESIII detectors, 
the detector response and digitization models, as well as the tracking of the detector running 
conditions and performances. An inclusive ip' MC sample is generated to study potential 
backgrounds. The production of the ip' resonance is simulated with the MC event generator 



KKMC 



15| . while the decays are generated with besevtgen [16|] for known decay modes 
with branching fractions being set at their PDG jsl world average values, and with lund- 



CHARM [17|] for the remaining unknown decays. The analysis is performed in the framework 
of the BESIII offline software system [l^ which provides the detector calibration, event 
reconstruction and data storage. 



III. EVENT SELECTION 

The selection criteria described below are similar to those used in previous BESIII anal- 
yses [qI, Q and are optimized according to the signal significance. 



A. Photon identification 

Electromagnetic showers are reconstructed by clustering EMC crystal energies. The 
energy deposited in nearby TOF counters is included to improve the reconstruction efficiency 
and the energy resolution. Shower identified as photon candidates must satisfy fiducial and 
shower-quality requirements. Photon candidates that are reconstructed from the barrel 
region (| cos^| < 0.8) must have a minimum energy of 25 MeV, while those in the endcaps 
(0.86 < I cos^l < 0.92) must have at least 50 MeV. Showers in the angular range between 
the barrel and endcap are poorly reconstructed and excluded from the analysis. To eliminate 
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showers caused by bremsstrahlung charged particles, a photon must be separated by at least 
10° from any charged track. EMC cluster timing requirements are used to suppress electronic 
noise and energy deposits from uncorrected events. The number of photon candidates 
is required to be 2 < A^^ < 10. 

B. Charged particle identification 

Charged tracks are reconstructed from hits in the MDC. For each track, the polar angle 
must satisfy | cos6\ < 0.93, and it must originate within ±10 cm from the interaction point 
in the beam direction and within ±1 cm of the beam line in the plane perpendicular to 
the beam. The number of charged tracks is required to be two with a net charge of zero. 
The time-of-fiight and energy loss dE/dx measurements are combined to calculate particle 
identification (PID) probabilities for pion, kaon, and proton/antiproton hypotheses, and each 
track is assigned a particle type corresponding to the hypothesis with the highest confidence 
level. Both charged tracks are required to be identified as kaons. 

C. Event selection criteria 

To choose the correct 77 combination for the vr^ or rj identification and to improve 
the overall mass resolution, a four-constraint kinematic fit (4C-fit) is applied under the 
hypothesis ip' — )■ '-yjK^K' constrained to the sum of the initial e^e~ beam four-momentum. 
For events with more than two photon candidates, the combination with the smallest is 
kept. Candidates with < 20 for this fit are retained for further analysis. Figure [T] shows 
the invariant mass distribution for the two selected photons. Signal candidates of vr" and r] 
mesons are clearly seen. 

1. Final selection of xp' — >■ K~^K~7V^ 

Candidates 7r° are selected by requiring the invariant mass of two photons, M^^, to satisfy 
the condition 0.117 GeV/c^ < M-y-y < 0.147 GeV/c^, an interval that is six times the vr" mass 
resolution (~5 MeV/c^). To suppress the background from ip' — )■ 'jXco, with Xco K^K~, 
it is required that the energy of the less energetic photon (E^^^^) is larger than 70 MeV. 
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FIG. 1: The invariant mass distribution for two photons in the selected ip' — >• 'yyK^K^ events. 

Background events from ip' — )■ -K^J/ip, with J/ip ^ K^K^, are removed by requiring that 
the mass of the two kaons satisfies \Mk+k- ~ > 7 MeV/c^, where mj/^ is the J /ip 

mass jst. 

There are in total 1158 ip' — > K^K^tt^ events selected from the data. A Dalitz plot of 
these events is shown in Fig. |2l Invariant mass spectra of tc^K^ and K^K^ are shown in 
Fig. [31 The two peaks in the Tt'^K^ mass spectrum correspond to the ir*(892)^ and K*j^, 
where K*j may be -ft'}(1430) or i^*(1410). A partial wave analysis (PWA), described below, 
is used to study the Dalitz plot structures. 

2. Final selection of ip' — >■ K~^K~r] 

The rj candidates are reconstructed using the two selected photons in '~f'jK^K~, and 
the T] yields are determined by a fit to the M^^ distribution. To suppress the background 
from ip' — 7- rjj/ip, with J/ip ^ K^K", the invariant mass of the two kaons is required 
to be less than 3.05 GeV/c^. The background from the decay ip' — )■ 7Xco/2, with Xco/2 
TT^ /rjK^K' {xci TT^K^K^ or rjK^K^ is forbidden), is suppressed by requiring that the 
lower energy photon should be outside of the range 115 MeV to 185 MeV. A Dalitz plot of the 
surviving events is shown in Fig. HI which is produced by using a loose rj mass requirement 
of 0.48 GeY/c < M^^ < 0.6 GeV/c^ compared to the mass resolution for r/ — )■ 77 (~7 
MeV/c^). The diagonal band shows a clean signal for ip' — )■ (pi] decays. 
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FIG. 2: The Dalitz plot for ip' K+R-tt^ 



IV. PARTIAL WAVE ANALYSIS OF tp' K+K'tt^ 



We perform a partial wave analysis of the decay tp' — )■ K^K 7r° in order to determine 
branching fractions for ip' K*(892)=^i^r'F and K*j^K^. 



A. The method 



The method of the PWA is similar to that utilized in a previous BES publication 20 |. 



The decay amplitudes are constructed using the relativistic covariant tensor amplitudes as 



described in Ref. 21|. For the decay ip' — ?■ K^K tt'^, the general form of amplitude reads: 

A{m) = ^P^{m)A^ = ^'^(m) ^ A,f/f , (3) 

i 

where %p^{m) is the polarization vector of ip' with a helicity value m; is the i-th partial- 
wave amplitude with the coupling strength determined by a complex parameter Aj. The 
differential cross section is given by 

^ = ^ E A^{rn)A*^{rn) = " (4) 

m=±l 
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FIG. 3: The invariant mass projection of the Dahtz plot (see Fig. [3|) for the xp' — )• K^K^n^ decay, 
(a) Mj^+^- is plotted with one entry per event, and (b) M^ox± is plotted with two entries per 
event. 

where P^- = P*^ = A^A*-?' and Fij = F*- = lY.l=iUtU*'' . Here, the sum over the -0' 
polarization is taken as m = ±1 since the ip' particle is produced from e~^e~ annihilation. 
The partial wave amplitudes Ui for the intermediate states, e.g. K*{892)^K^, K;{U30)^K^ 
etc., are constructed from the K~^, K~ and vr'' four- momenta. In the amplitude, the line 
shape for the resonance is described with a Breit-Wigner function: 

1 



BW{s) 



(5) 



M^-s- iMT ' 

where s is the invariant-mass squared, and M and F represent the mass and width, respec- 
tively. 

The relative magnitudes and phases for amplitudes Ui are determined by an unbinned 
maximum likelihood fit. The joint probability density for observing the events in the data 
sample is 

N 

£ = l[P{x,), (6) 



i=l 
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FIG. 4: The Dalitz plot for -0' rjK+K^. 

where P( probability to produce event i with four-vector momentum Xi 

{PK+,PK-,PTTo)i- The normahzed P{xi) is calculated from the differential cross section 

{da/d^)i 



PiXi) 



(7) 



where the normalization factor aMC is calculated from a MC sample with Nmc accepted 
events, which are generated with a phase space model and then subject to the detector 
simulation, and are passed through the same event selection criteria as applied to the data 
analysis. With an MC sample of sufficiently large size, the (Tmc is evaluated with 



MC 



i=l 



For technical reasons, rather than maximizing £, S* = — In £ is minimized using the package 
FUMILI 



B. Background subtraction 

The number of non-7r° background events in the selected K^K^tt^ data sample, estimated 
from a 7r° sideband defined by M^^ e [0.079,0.109] and [0.165,0.195] GeV/c^, is 43±7 
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events. The MC simulation shows that these background events are mainly due to if)' — )■ 
IXcJi Xcj 'jK^K^ or n^K^K^ . A low level of non-K~^K^ background (3 events) comes 
from ip' — > TT^TT^J/ip, J/ip ^ fi'^fi' due to a misidentification of muons as kaons. 

Events from the QED process, e+e^ — > 7* — i- K^K^ir^ produced at a center-of-mass 
energy corresponding to the mass of the ip' peak, have the same final state as our signals of 
interest. Background from this source is estimated from two data sets taken at y/s = 3.773 



0, 



the events 



GeV and 3.65 GeV. Since the decay of '?/'(3770) — )■ K^K^n'^ is not observed 
obtained at ^/s = 3.773 GeV are regarded as all due to the QED process. After normalizing 
their integrated luminosities to that of the ip' sample, the number of events obtained at each 
of the data sets are 195±3 and 195±27, respectively, and in good agreement with each other. 

The QED background events at the ip' peak are generated using a model determined by 
performing a PWA fit to the data set taken at 3.773 GeV. As a cross check, the model with 
the determined coupling strengths is used to generate MC samples and compared with the 
data set taken at 3.650 GeV. Figure [5] compares mass distributions obtained from MC events 
with those obtained from experimental data. Here MC and experimental data were generated 
or taken at ^/s = 3.650 GeV. For the K^K~ and Kt^^ invariant mass distributions, the data 
and MC agree well within statistical errors, and a peak around Mt^qk± = 1-4 GeV/c^ can 
be seen. 

In the PWA fit, background events obtained from MC simulation or vr" mass sideband 
are used to account for the background events in the data using a negative log-likelihood 
value. Hence, the complete log-likelihood function is: 

In £ = In Cdt - X] In Cbg, (9) 

where Cat and C^g are the likelihoods determined with the data and background events, 
respectively. The backgrounds are divided into two kinds: reducible background and irre- 
ducible background (QED background). This technique of background treatment assumes 
no interference between signal and irreducible background events. This method has been 



used in the analysis of Crystal Barrel data 23|] and BESII data 20|, |24 |. 
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FIG. 5: The K^K^ (one entry per event) and Kit (two entries per event) invariant mass distri- 
butions at ^/s = 3.65 GeV. The dots with error bars are data and the histograms are MC events 
as described in the text. 

C. Analysis results 



Motivated by the structures seen in the Dalitz plot (Fig. [2]) and its projections (Fig. |3]), 
the decay modes listed in Tables HI and HTl are considered in the PWA fit. Only the modes 
with a statistical significance larger than 5 standard deviation (a) are taken as the best 
solution, which includes the resonances i\:*(892)±, ir*(1430)±, i\:*(1680)± and p(1700), and 
the non- resonance mode K^K^tt^ (see Table [I]). The significance of a mode is calculated 
by comparing the difference of the S{= — In C) values between the fit with and without 
that mode. The non-resonance mode is described as a P— wave K^K' system. For the 
charge-conjugate channels, the coupling strengths in amplitudes are the same. Each mode 
in the amplitude introduces two parameters are determined by the PWA fit, the magnitude 
of the coupling strength and the phase angle. 

Other intermediate states, like p(770), p(1450), p(1900), p(2150) in the K^K^ final states, 
and ii'*(1410)^ and i^'*(1980)^ in the i^^K^ final states, were considered and tested in the 
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PWA fit. Adding them to the best solution does improve the fit quality, but these additional 
modes have a statistical significance of less than 5a (see Table HTj) . The p(770) may decay to 
K^K~ if its mass is larger than K^K^ threshold, but its significance is 4. 60". A P— wave 
Tt^K^ system as an additional non-resonance contribution was tried and had a significance 
of 1.9cr. The variations to the iC*(892)^ and -ft'2(1430)^ signal yields by including these 
intermediate states are included as a systematic uncertainty. 

TABLE I: The significance and number of events of each resonance under the best solution. 



Decay Fitted events Significance(o") 
K*{m2)^K^ 224±21 26.5 
K*{IA2,Q)^K^ 251±22 21.0 
i^*(1680)±E:=F 115±20 11.1 
p°(1700)7r° 59±10 8.7 
K+K--K^ 721±60 18.8 



TABLE II: Significance for additional resonance. 



Decay Sij 


;nificance(o") 


pO(770)7rO 


4.63 


pO(1450)^o 


4.40 


p°(1900)7r0 


1.13 


p0(2150)7r0 


3.21 


p!](1690)7r° 


1.84 




2.23 




2.14 




3.05 




3.26 


non-resonance {K^tt^) 


1.89 



For intermediate states around Ktx invariant mass of 1.43 GeV, there are four established 
resonances, namely, ifi(1400), if*(1410), i^o(1430) and i^2(1430); according to the spin- 
parity conservation, only ^^2(1430) and K*{1A1Q) are allowed. If (1430)^i^'^, which is 



15 



the best solution in the PWA, is replaced with i^*(1410)^-ft'^, the fit fails to match the 
data, and the log-likelihood gets worse by 126, and the contribution from the ir*(1410) is 
negligible. If K*{U10)^K^ is taken in addition to K^{U30)^K^ to the best solution, the 
log-likelihood only improves by 3.65, corresponding to a significance of 2.2a. 

The non-resonance decay ip' — K^K^tt^ is indispensable in the fit, with a statistical 
significance of 19a. We have tried to replace it with a broad resonance, such as p(2150)7r°. 
The fit fails to match the data, and the log-likelihood gets worse by 95. Note that the total 
number of fitted events 1370 ±70 in Table [I] is larger than the number of net K^K^tt^ events 
917(=1158-241) due to the destructive interference among the included resonances. 

The numbers of fitted events given in Table [T] are derived from numerical integration of 
the resultant amplitudes as done in Ref. {2^. The statistical errors are derived from the S 
distribution versus the number of fitted events; one standard deviation corresponds to the 
interval that produce a change of log-likelihood of 0.5. When performing the PWA fit to 
the data, the masses and widths of the intermediate states are fixed at the PDG values, and 
their errors quoted in the PDG are used to estimate the associated systematic errors. 

Figure depicts a comparison between the data and the best solution obtained from the 
PWA fit to the data. Here the projected M^+k- and M^o^± mass distributions are shown. 
They are in general in a good agreement except for several points at the low M^+k- mass 
region. An additional p(1450)'7r'' to the best solution in the PWA helps to improve the fit 
quality through destructive interference (see Fig. [7]). The statistical significance of this 
additional mode is only about 3.2cr and it only brings a small difference in signal yields, 
3.3% for K*{892)^K^ and 0.4% for K2{U30y^K^. These yield differences are taken as a 
systematic uncertainties to account for additional resonance contributions to the low Mj^+k- 
mass region. 

The goodness of the global fit is determined by calculating a xin defined by 

xln = E Xl with x| = E ^ ^ , (10) 

j=i 1=1 i« 

where N^'^ and Nf/^ are the number of events in the i-th bin for the distribution of the ?'-th 
kinematic variable. If the measured values ^fi^ are sufficiently large, then the Xaii statistic 
follows the distribution function with the number of degrees of freedom (ndf) equal to 
the total bins of histograms 27| minus the number of fitted parameters; and the individual 
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FIG. 6: The results of fit to (a) K^K^ and (b) t:^K^ mass distributions for the data, where points 
with error bars are data and histograms are total fit results. The dashed histograms are the sum 
of the background sources, including QED and non-X"'"i^~7r'^ contributions. 

x| gives a qualitative measure of the goodness of the fit for each kinematic variable. 

For the 3-body decay ip' — )■ K^K^tt^, there are 5-independent variables, which are se- 
lected as the mass of the K^K~ system {Mx+k-)i the mass of the t^^K^ system (M^ox±), 
the polar angle for the vr*^ (^jro), the polar angle for the K~ (Ok-), and the azimuthal angle 
for the {(f)x+), where the angles are defined in the tp' rest frame. Figure [8] compares the 
angular distributions between the best fit solution and the data, and a good agreement can 
be observed. A sum of all these Xj values gives xlw = 147.70, and the total number of de- 
grees of freedom (126) is taken as the sum of the total number of bins having non-zero events 
minus the total number of parameters in the PWA fit. The global fit goodness xln/^df is 
1.2. 
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FIG. 7: The results of fit to (a) K^K~ and (b) tt^K^ mass distributions for the data; where points 
with error bars are data; histograms denote total fit results with an additional mode of /9(1450)7r'' 
being added to the best solution of the PWA fit (see Fig. [6]) . The dashed histograms are the sum 
of the background sources, including QED and non-X"'"i^~7r^ contributions. 

D. Branching fractions 



Branching fractions for ^' ir*(892)+ir" + c.c, f K^{1A30)+K- + c.c, and the 
inclusive decay ip' — )■ K^K^tt^ (including all resonances) are calculated 



Br{ilj' K*+K- + c.c.) 



eN^>Br{K*+ ir+7rO)fir(7rO ^ 77)^ 

TVfobs 



(11) 



eN^'Br{n^ 77) ' 

Here Br{K*^ K^t:"^) is the branching fraction for ir*(892)+ (33.23%) or K*(1430)+ 
(16.60%) resonances; A'^^^''^ is the signal yield obtained from the PWA fit (224 ± 21 and 
251 ± 22 for ir*(892) and ir2*(1430), respectively); iV^+^^-^o is the net number of K+R-n^ 
events (917±37); N^i = (106 ±4) x 10^ is the number of ip' events[25!]; and e is the detection 
efficiency. To determine e, the intensity from the amplitudes is used to weight both the 
complete set of generated MC events and the set which survives the selection procedure, 
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FIG. 8: The fit results of the angular distributions, where the points with error bars are data and 
histograms are the fit results, (a) cos^? distribution for vr^, (b) cos^? distribution for K~^, and (c) 
sin^ distribution for K~^, here the angles are defined in the tp' rest frame. 



and the ratio between these two weighted sets is taken as the detection efficiency. 
The branching fractions are measured to be: 

Br{ip' TT^K+K-) = (4.07 ± 0.16) x 10'^ (12) 
Briijj' K*(892)+K- + c.c.) = (3.18 x 0.30) x 10'^, (13) 
Briil^' K;{U30yK- + c.c.) = (7.12 ± 0.62) x 10-^ (14) 



where the errors are only statistical. 
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FIG. 9: The K^K^ invariant mass selected in ip' — t- ^^K^ . 
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FIG. 10: (Color online) The invariant mass distribution of two photons in the selected ip' — )■ 
077 events; the solid line shows a fit to vr*'; the dashed line shows the fitted background and 
comparison to the backgrounds estimated with (j) sideband (line histogram) and MC simulation 
(dashed histogram). 

V. V' 7r°0 

The (f) candidates for tp' — (pii^ are reconstructed using the two kaons selected in the 
decay ip' — )■ K'^ K" '-yy . Figure |9] shows the invariant mass distribution of the two kaons, and 
a signal is clearly seen. The candidates are selected by requiring \Mk+k- — < 10 
MeV/c^, where M^+k- and are the invariant mass of the two kaons and the mass of the 
jsi- Background sources from the initial state radiation process e+e~ — )■ 70 are suppressed 
by requiring that the energy for the energetic photon is less than 1.6 GeV. Figure [TO] shows 
the invariant mass distribution of the two photons after the selection criterion is applied. 
No significant 7r° signal is observed. 

The number of observed events for if)' — ix^cj) is obtained by fitting the mass distribution of 
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the two photons as shown in Fig. [TOl The hne shape of vr*^ is taken from the MC simulation, 
and the background shape is taken as a first-order Chebychev polynomial function. The fit 
results are shown in Fig. [10] and the significance of vr'' signal is less than 3. Oct. The upper 
limit of observed events is estimated using the Bayesian approach to be N^^ = 6 at the 
90% confidence level. 

The upper limit on the branching fraction for i/)' — )■ 7r^(f) is calculated with 

Bri^P' -> vr 0) < ^jy^,5^(^o ^ ^ K+K-){1 - a^v^) ' ^^^^ 

where Br(j{^ — )■ 77) and Br{(j) — )■ K^K^) are the branching fractions for tx^ — )■ 77 and 
(f) — 7- K~^K^ , respectively; N^/ = (106 ± 4) x 10^ is the number of total tp' decays; e = 
35.63% is the detection efficiency that was determined using MC events generated with the 
angular distribution 1 + cos^ 6 for ip' — )■ 7r°0, where 6 is the (p polar angle, a^'^'^ = 5.8% 
is the systematic error as listed in Table IIIII The upper limit of the branching fraction is 
Br{tp' (pn°) < 4.0 x 10"^ at the 90% C.L. 

VI. ip' rjK+K- 

A. Background analysis 

Background sources for ip' — )■ rjK^K^ are studied with the ip' inclusive MC sample. The 
dominant background comes from ip' — )■ •y^psnK^K' , where 'Jfsr is a final-state radiation 
photon, Ip' — )■ 7Xc2; with Xc2 K^K^tx^ and K^K'rj. The MC simulation shows that the 
M^^ mass distribution of sum of these events in the region of the t] meson is a smooth and 
well modeled with a polynomial function. 

Background events from QED processes are studied using events taken at ^/s = 3.773 
GeV that are selected with the same criteria applied to the ip' data. The signal yields 
are extracted with the same fit procedure used for the ip' data. For r/^, the contribution 
from the resonance decay ip{3770) — )■ r](p is estimated to be 450 ± 112 events using the 



measured cross section a = 2.4 ± 0.6 pb 26|]. After subtracting the resonance decays, the 



QED yield for the e+e^ — t- r](f) at y/s = 3.773 GeV is determined to be 268±115 events. For 
TjK^K^, the observed events are considered to be exclusively from QED processes because 
the ip{3770) — t- rjK^K^ has not observed jst. At the ip' peak, the QED background sources 
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are estimated to be 16±7 events for the T](f) and 4±1 events for the riK^K~ according to 
the luminosity normahzation. As a cross check, we use the data taken at ^/s = 3.65 GeV to 
determine a QED background of 25 ± 9 events. The difference between the two estimates is 
taken as a background uncertainty and included into systematic errors. 

B. Fit results 

We performed a two-dimensional unbinned fit to the scatter plot of M^+k- versus M^^ 
distribution assuming that M^+k- and M^^ are independent variables. Motivated by the 
structures seen in the Mk+k- distribution, resonances including 0(1020), 03(1850) and 
0(2170) are added to the fit. The fit function includes the line shapes describing the two-body 
decays ?70(1O2O), ?703(185O), ?70(217O), the non-resonant decay r]K~^K~ , and the background. 
The f] line shape is obtained from a MC simulation; the line shapes for the 0(1020), 03(1850) 
and 0(2170) are described as non-relativistic Breit-Wigner functions with their masses and 
widths fixed to the PDG values. The Breit-Wigner function of all the states are convolved 
with a detector resolution function. The background shapes for the M^^ and the M^+k- 
mass distributions are taken as first- and second-order polynomials, respectively. 

The fit results after projecting to the mass distributions are shown in Figs. [TT] and O 
The signal yield for the tp' — )■ 7^0 channel is 232 ± 16 events. Adding the 03(1850) and 
0(2170) resonances to the fit improves the fit quality with a statistical significance of 3.8a 
for the 03(1850), and 3.1a for the 0(2170). The goodness of the fit is x^/ndf = 0.32(0.43) for 
the M^^{Mk+k-) distribution. The yields of ?703(185O) and 7^0(2170) plus the contribution 
from the non-resonance decay ip' — ?■ rjK^K^ totals 288 ± 27 events. After subtracting the 
QED background, the net signals are 216 ± 16 events for ip' — )■ 7^0, and 284 ± 27 events for 
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FIG. 11: (Color online) Fit results projected to the two-photon invariant mass distribution M^^. 
Dots with error bars are data. The solid line is the total fit results, and the dashed-dotted and 
long-dashed lines are the results of ijcj) and rjKK contributions, respectively. The short-dashed line 
is the background contribution. 

C. Branching fractions 

Branching fractions are calculated from the relations 



Here N^'^j^ = 284±27 and A^^^* = 216±16 are the numbers of net signal events; Br{ri — )■ 77) 
and Br{(f) — )■ K^K~) are the branching fractions for the — )■ 77 and — )■ K^K~ decays, 
respectively; er,KK = 22.10% and Er,^ = 33.53% are the detection efficiencies determined from 
MC simulations, whose angular distributions match the data; ErfKK is a weighted average 
for ip' — riK'^K~, 7703(1850) and 770(2170). The branching fractions are calculated to be 
Br{ij' r]KK) = (2.97 ± 0.28) x 10"^ and Br{f r]4>) = (3.08 ± 0.29) x IQ-^, where the 
errors are only statistical. 

VII. SYSTEMATIC ERRORS 

The systematic errors in the branching fraction measurement originated from following 
sources are considered: 




(16) 



Br{ilj' — > 7/0) 




(17) 
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FIG. 12: (Color online) Fit results projected to the K^K~ invariant mass distribution Mj^+j^- 
for (a) the 0(1020) resonance, (b) the (/)(1850) and (/'(2170) resonances. Dots with error bars are 
data. The solid lines are the total fit results, and the dashed-dotted and long-dashed lines are the 
results of ri(f) and ijK^K^ , respectively. The short-dashed line is the background contribution. 

1. photon efficiency 

The soft and hard photon efficiencies are studied using tp' — ti^ti^ J /ip, J /ip — )■ 
e~^e~,^~^^~ and J/ip ^ pn ^ tt^ti^tt'^ decays. The difference in the photon effi- 
ciency between the MC simulation and data is 1%, which is taken as a systematic 
uncertainty. 

2. kaon tracking and PID efficiency 

The uncertainties of kaon tracking and PID efficiency are studied using a sample of 
J/ip ii'*(892)°K^ + c.c. K^K^TT' + c.c. A'+tt^tt+tt" + c.c. events as done in 



19j. The uncertainties for both tracking and PID are determined to be 1% per track. 



3. Number of ip' events 

The number of ip' events is determined using its hadronic decays. The uncertainty is 
4% [251. 



4. branching fractions 

The uncertainties of branching fractions for i^'*(892)^/i^r2(1430)^ — > K^tc^, n^/i] — )■ 
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77 and — i- K^K are taken from the world average values [Sj. 

5. kinematic fit 

The differences between the MC simulation and data in the distribution of the 
kinematic fit arise mainly due to inconsistences in the charged track parameters. The 
kaon track parameters in the MC simulation are corrected by smearing them to match 
the data. The difference in the detection efficiency between with and without making 
a correction to the MC is taken as a systematic error. The uncertainties are listed in 
Table [ml 

6. the 7r° mass window 

The uncertainty due to the 7r° mass window is studied by comparing the 7r° selection 
efficiency obtained in the MC and the data. The uncertainty is 1.1%. 

7. fit uncertainty 

The fit uncertainties in the rjK'^K' and r^^ modes are determined by changing the 
fit range and background shapes. The fit range of two photons is changed to be [460, 
620] MeV/c2 or [470, 670] MeV/c^. It is estimated to be 3.6% (0.6%) for riK+R- 
(?70). The background function is changed from Ist-order to 3rd-order polynomials. 
The uncertainties due to the background shapes are 1.6% and 0.4% for rjK^K" and 
rjcp, respectively. 

8. QED backgrounds 

The QED background subtracted from 7^0 is determined with the data taken at ^/s = 
3.773 GeV and at ^/s = 3.65 GeV. The difference in the number of QED events 
between these two samples is 4.5%, which is taken as the QED background associated 
uncertainty. 

9. additional resonances for rjK^K^ 

The existence of 03(1850) and 0(2170) intermediate states in riK^K~ cannot be de- 
termined due to the low statistics. The difference between the branching fractions 
determined by including and excluding these two resonances is taken as a systematic 
error of 4.0%. 

All above systematic errors are listed in Table IIIII 
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For K^K 7r°, the uncertainties from the PWA fit are hsted below: 



1. Breit-Wigner form 

The uncertainty due to the resonance hne shape is evaluated by using the Breit-Wigner 
function with a width r(s) dependent on the energy, i.e. 

p{s) 



s 



2 r l2L+l 



(18) 



where s is the resonance mass squared; m and Fq are the nominal mass and width, 
respectively; p{s) is the magnitude of resonance momentum; L is the angular mo- 
mentum for the ip' decays into a two-body final state. The differences between the fit 
yields determined with a constant and an energy-dependent width are taken as system- 
atic errors. They are evaluated to be 0.1% and 0.9% for the A"* (892) and A'2*(1430), 
respectively. 

2. additional resonances 

The uncertainties from additional resonances, listed in Table [Tll are determined by 
adding them to the best solution of PWA fit one-by-one. The differences between the 
fit yields determined with and without the additional resonance are taken as systematic 
errors. For the non-resonant mode ip' — K^K^tt^, the uncertainty due to the P-wave 
K^K~ system in the PWA fit is evaluated by replacing it with a P-wave Kt{ system. 
The difference in the fit yields is taken as a systematic error. 

3. Yio\i-K^ K~Ti^ background 

The number of won-K^ K~ti background events is obtained from a vr^-sideband analysis 
and an exclusive MC simulation. The difference in the signal yields corresponding to 
one standard deviation of this background is taken as a systematic error. 

4. the QED background 

The QED background used at y/s = 3.686 GeV is produced via a MC simulation with 
amplitude information obtained from a PWA fit to the data taken at ^/s = 3.773 GeV. 
The uncertainty is estimated by replacing this QED background with the continuum 
data taken at y/s = 3.65 GeV. The difference of the fitted yields between these two 
approaches are 0.8% and 9.9% for ^^2(1430) and -A* (892), respectively, and used as 
systematic uncertainties. 
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5. uncertainty of -ft'*(1680) and p(1700) widths 

The decay widths of -ft'*(1680) and p(1700) have large uncertainties; the world average 
values are r;^.(i68o) = 322 ± 110 MeV and rp(i7oo) = 250 ± 100 MeV fl- The sig- 
nal yields were re-obtained using widths that are changed by one standard deviation 
with respect to the nominal value. The differences in signal yields between these two 
methods are taken as systematic errors. 

6. uncertainties of masses and widths for the K*{892) and K*{1A30) 

In the PWA fit, the masses and widths for the K*{892) and K*{U30) are fixed to 
the world average values. The differences in fit yields obtained by changing these 
parameters by one standard deviation are taken as systematic errors. 

All systematic errors from the PWA fit are listed in Table IIVI 

Combining the systematic uncertainties from the PWA fit and the n^K~^K~ event se- 
lection gives total systematic errors of ^1'^% and ^8^^% for ip' — )■ K* {892)^ K~ + c.c. and 
K;{U30yK- + c.c, respectively. 



VIII. SUMMARY AND DISCUSSION 



Using (106±4) x 10^ ip' decays accumulated with BESIII, we measured branching fractions 
for the ijj' K*{892)+K- + c.c, K*{U3f))+K- + c.c, r](j), TT°(f),TT°K+K- , and 7]K+K- 
decays. The helicity forbidden decay ip' — > (1430)+/^^ + c.c. is observed for the first 
time, and its branching fraction is measured; this reflects a violation of the helicity selection 
rule m. Table |V] gives an overview of our results with comparisons with BESII- and CLEO- 
measurements and world average values. The precision of our measurements is better for 
all the modes, including a tightened upper limit for 7r°0. In the measurement of Br{ilj' — )■ 
n^K'^ K~), all intermediate states are included in the branching fraction, while for the 
measurement of Br{il)' — K~^K~ri), ip' — > r/^ is excluded. The measurements of branching 
fractions for the ip' — )■ K*{892)^K^ + c.c. and ?70 are consistent with BESII results within 
la, and CLEO measurements within 2a. 

Using the world average values of branching fractions for J/ip decays, the Qh values 
are calculated and listed in Table |Vl For ip' — )• K*{892)^K^ + c.c and rjcp, the Qh values 
significantly deviate from the expected value of 12%. 
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TABLE III: Summary of all systematic errors (%). 
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2 


PID 
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2 


Kinematic fitting 


1.9 


3.2 
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Number of tjj' decays 
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5.8 



TABLE IV: Summar>- of s>-st (^nia1ic uncc^rtaintic^s from tlic^ VWA {%). 



Sources 


K*(892)±if=F K*{Um)^K^ 


Breit-Wigncr 


-0.1 




+0.9 


Additional states 


+5.2 

-6.9 




+10.3 
-4.6 


Non-if+K~"7r background 


+1.4 
-1.6 




+1.2 
-1.0 


QED background 


-0.8 




+9.9 


K*(1680), p(1700) width 


+0.5 
-1.1 






-2.0 


K*(892), ii:2*(1430) Mass and width 


+0.4 
-0.4 




+0.3 



Total 


+5.4 
-7.3 




+14.3 
-5.1 
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n 

TABLE V: Summary of the measured branching fractions compared with PDG [5] values, together 
with CLEO 3| and BESII [4,] measurements. The upper Umit is given at the 90% confidence level. 
The first error is statistical, and the second error is systematic. Here e, N"'''', and Br denote the 
detection efficiency, the number of observed events, and the branching fraction, respectively. The 
variable Qh is defined by Eq. ([1]). 



Mode(?A' ^) 


e(%) N"*-^ Br(xl0-5) Q,, (%) 


PDG(xl0-'^) CLEO(xlO-5) BESII(xlO-^) 


tt'-'K^K^ (inclusive) 


21.52 917±37 4.07±0.16 ± 0.26 — 


< 8.9 [2j 




K* (892)+ K- +C.C. 


20.25 224±21 3.18±0.30lJ^j^ 0.62 ±0.09 


1 7+0.8 
^■'-0.9 


1.3 ± 1.0 ±0.3 2.9 ± 1.3 ±0.4 


KI{U30)+K- +C.C. 


20.28 251±22 7.12±0.62l J'^J > 2 






T]K^K^ {T](j) excluded) 


22.10 284±27 3.08±0.29 ± 0.25 — 


< 13 


< 13 — 




33.53 216±16 3.14±0.23 ± 0.23 4.19 ±0.61 


^•^-0.8 


2.0± 1.1±0.4 3.3± 1.1±0.5 




35.63 <10 <0.04 - 


< 0.4 


0.7 0.4 
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